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ABSTRACT 

The mass of the Galactic dark matter halo is under vivid discussion. A recent study by Xue et al. (2008, 
ApJ, 684, 1 143) revised the Galactic halo mass downward by a factor of ~2 relative to previous work, based 
on the line-of-sight velocity distribution of ~2400 blue horizontal-branch (BHB) halo stars. The observations 
were interpreted in a statistical approach using cosmological galaxy formation simulations, as only four of the 
6D phase-space coordinates were determined. Here we concentrate on a close investigation of the stars with 
highest negative radial velocity from that sample. For one star, SDSSJ153935.67+023909.8 (J1539+0239 for 
short), we succeed in measuring a significant proper motion, i.e. full phase-space information is obtained. 
We confirm the star to be a Population II BHB star from an independent quantitative analysis of the SDSS 
spectrum - providing the first NLTE study of any halo BHB star - and reconstruct its 3D trajectory in the 
Galactic potential. J1539-I-0239 turns out as the fastest halo star known to date, with a Galactic rest-frame 
velocity of 694;^22? ^^ ^ ' (^^^^ uncertainty range from Monte Carlo error propagation) at its current position. 



The extreme kinematics of the star allows a significant lower limit to be put on the halo mass in order to keep 
it bound, of Mhaio ^ ^-"^-I'l ^ 10'^ Mq. We conclude that the Xue et al. results tend to underestimate the true 
halo mass as their most likely mass value is consistent with our analysis only at a level of 4%. However, our 
result confirms other studies that make use of the full phase-space information. 



Subject headings: dark matter — Galaxy: halo — stars: atmospheres 
kinematics and dynamics — stars: Population II 
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1. EVITRODUCTION 

Knowledge of the properties of dark matter halos is an im- 
portant issue for our understanding of galaxy formation and 
evolution, and for unveiling the nature of dark matter. The 
halo of the Milky Way therefore is of highest interest, as 
it allows unique observational c onstraints to be obta ined for 
testing theoretical models (e.g. iNavarro et al.l 119961) . Sev- 
eral observational campaigns - e.g. the Sloan Digital Sky 
Survey (SDSS , York et al. 2000), the RAdial Velocity Ex- 
periment (RAVE, ISteinmetz et alJ i2006), the Sloan Exten- 
sion for Galactic Understanding and Exploration {SEGUE, 
lYannv et"al] l2009h - provide the tracers for studying halo 
properties, like the total mass of the halo and its extent. 

Several studies in the past decade determined the halo mass 
from ever increasing samples of halo stars, globular clusters 
and/or satellite galaxies. However, it is in fact only a few 
objects at highest velocity that are affecting largely the mass 
estimates (.Sakamoto et all 120031: ISmith et al.ll2007h . While 
larg er halo masses of about 2 x lO'^Mp were favou red ear- 
lier dffilkinson & Evans"1999l IS^amotp et al."2003'), lower 
masses of a bout half this val ue were d erived more recently 
(iBattagha et alJ2005l:ISmith et al. 2007; Xue et al. 2008). The 
precise value determines, e.g. among others, whet her satellite 
galax i es like the Magellanic Clouds (e.g. Kal liyavalil et alj 
l2006t ICosta et al] 1200^) or hyper- velocity stafl candidates 
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' Hyper-velocity stars move at such high velocity that they may be gravita- 
tionally unbound to the Galaxy. Their supposed place of origin is the Galac- 
tic center, where they may have been accelerated by gravitational interactions 



dAbadi et al.ll2009h are on bound orbits, or not. 

Most of the previous studies had to rely substantially on 
the distributions of radial velocities in their samples to de- 
rive their conclusions as full space motions for halo objects 
are unavailable in many cases at present. In such cases only 
four coordinates (i.e. two position values, distance and ra- 
dial velocity, RV) of the 6D phase space are determined and 
the missing data (the proper motion components) are handled 
in a statistical approach. E.g., radial velocities for more than 
10,000 blue halo stars from the SDSS were m easured in the 
yet most extensive study by IXue et all (120081) . The sample 
was composed of blue horizontal-branch, blue straggler, and 
main-sequence stars with effective temperatures roughly be- 
tween 7,000 and 1 0,000 K according to their colours. Here, 
we focus on the fastest stars of the Xue et al. sample in terms 
of negative line-of-sight velocity, indicating a bound orbit. 

For one of them we were able to measure a significant 
proper motion, which allowed a detailed 3D kinematic in- 
vestigation when combined with a quantitative spectroscopic 
analysis that facilitated the determination of the star's dis- 
tance. SDSSJ153935.67+023909.8 (J1539-(-0239 for short) 
is an inbouncQ Population II horizontal branch star with a 

with the super-massive black hole (Hills 1988). 

^ An extragalactic origin of the star as an unbound low-mass hyper- 
velocity star from another galaxy is imaginable but unlikely. The local vol- 
ume is devoid of galaxies in a wide range (several lOdeg) around the infall 
direct ion of J1539+0239 at present (for the Local Group see van den BergH 
119991 and later discoveries). If ejected early after its formation, the star could 
have travelled ~8 Mpc during its lifetime, such that a final answer on this may 
be only obtained when full space motions of the neai'by galaxies are known. 
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Galactic rest-frame (GRF) velocity of ~700kms"' at its cur- 
rent position, making it the fastest halo object known. This 
allows a significant lower limit for the total halo mass of the 
Galaxy to be set. The present work provides a glimpse to the 
detailed kinematic investigations feasible once the European 
Space Agency's Gaia satellite mission (e.g. lTuron et al.l2005l) 
becomes operational, at much higher precision. 

2. TARGET SELECTION AND PROPER MOTION 



In a previous paper ("T illich et al. r2Q09') we already inves- 
tigated the high-velocity tail of the Xue et al. (2008) sam- 
ple and found a hyper-velocity candidate of spectral type 
A among the stars with highest positive radial velocities in 
the GRF. Here we study the most extreme stars approach- 
ing us, applying the same techniques. We selected all stars 
with GRF velocities vqrf < -350 km s ' from the RV-based 
sample of Xue et al.l ('2008) and obtained 5 targets for which 
we attempted to measure proper motions. All available in- 
d ependent positio n measuremen ts on Schmidt plates (APM 

- IMcMahon et aH dlOOa) : SSS - IHamblv et al.1 (llOOl ) were 
collected and co mbined with the SPSS and other ava ilable po- 
s itions(CMC 14lCarlsberg-Meridian-Catalog'('2006'); 2MASS 

- ICutriet al.l (12003); UKIDSS - Lawrence et al. (2007)) for 
a first linear proper motion fit. However, there were even 
more measurements of Schmidt plates, from up to 14 dif- 
ferent epochs in case of overlapping plates of the Digitised 
Sky Surveyfl FITS images of 15 by 15 arcmin size were 
extracted from all available plates and ESO MIDAS tools 
were used to measure positions. For this purpose, we se- 
lected compact background galaxies around each target, iden- 
tified from SDSS, to transform the target positions on all the 
Schmidt plates to the SDSS system. The small fields allowed 
us to apply a simple model (shift-i-rotation) and to achieve 
a higher accuracy in our final proper motion fit for all our 
targets (see e.g. Fig. [T]i. We consider a positive detection 
if the proper motion errors cr of both components are be- 
low 3.0 mas yr ' and if at least one of the components is sig- 
nificant (above 3cr). For J1539H-0239, the brightest of our 
five targets, we found a highly significant proper motion of 
/i„cos 6 - -10. 6+1. 6 mas yr"' and/jg = -10.0+2.3 mas yr"', 
whereas for the other four fainter stars (with V magnitudes 
between 17 and 20) the proper motion was consistent with 
zero. The proper motions as well as their transformations 
from equatorial to Galactic coordinates are summarized in Ta- 
ble H 

3. OBSERVATIONS AND QUANTITATIVE SPECTROSCOPY 

In order to exclude RV variability, we re-observed 
J1539+0239 with the TWIN spectrograph at the 3.5m tele- 
scope on Calar Alto in May 2009. Radial velocities were 
derived by ;k'^ -fitting of adequate synthetic spectra over the 
full spectral range, yielding a heliocentric radial velocity of 
Vrad - -366.6 ± 4.0km s"' for the TWIN spectrum, which is 
consistent with the Vrad - -372.6 + 5.8 km s"' measured from 
the SDSS data within the mutual uncertainties. We use the 
latter value for the kinematic study. 

A quantitative analysis of J1539+0239 was carried 
out following the hybrid NLTE approach discussed by 
[Frzvbilla et al. (2006). In brief, line-blank eted LTE mode l at- 
mospheres were computed with ATLAS9 dKuruczil 1993b and 
NLTE (and LTE) line-formation calculations were performed 
using updated versions of Detail and Surface dGiddingsl 

^ http://archive.stsci.edu/cgi-bin/dss_plate_finder 
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Fig. 1. — Linear fit of the position measurements for J1539+0239 with a 
zero epoch of 1975.39. 



198iVButle r & Giddin gs 1985). Many astrophysically impor- 
tant chemical species were treated in NLTE, using state-of- 
the-art model atoms (H: Przybilla & Butler 2004a; C i: Przy- 
billa et al. 2001b; Ni: Przybilla & Butler 2001; Oi: Przy- 
billa et al. 2000; Mg i/ii: Przybilla et al. 2001a; Ti ii and Fe ii: 
Becker 1998). 

The effective temperature Teff and the surface gravity log g 
were determined by fits to the Stark-broadened Balmer and 
Paschen lines and an ionization equilibrium, here of Mg i/ii, 
in analogy to previous work on hyper- velocity sta rs at similar 
temperatures (iPrzvbilla et d]|200l IMich et al.ll2009h . The 
stellar metallicity was derived by model fits to the observed 
metal line spectra. Results are listed in Table [T] and a com- 
parison of the resulting final synthetic spectrum with obser- 
vation is shown in Fig. |2] Overall, excellent agreement is 
obtained for the strategic spectral lines throughout the entire 
wavelength range. Our stellar parameters (reff= 7700+250 K, 
log^ = 3.00 + 0. 15) are con sistent with those derived in the 
LTE analysis bv iXue et al.1 (12008) : Teff^ 7807 K, \ogg = 
3.16. We constrained the errors in the stellar parameters by 
the quality of the match of the spectral indicators within the 
given S/N limitations. 

Its parameters place J1539+0239 on the horizontal branch 
at a mass of 0.68+0.05 Mq as derived by comparing the posi- 
tion of the star in the T^ff -log ^-diagram to pr edictions of the 
evolutionary models of Dorman et al.l (Il993h . No rotational 
broadening was detected at the resolution of the SDSS spec- 
trum. The metallicity is lower than solar by a factor of 100 
and the abundances of the a-elements are enhanced by about 
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Fig. 2. — Compaiison of NLTE spectrum synthesis (thick line) with observation {SDSS, thin wiggly line) for J1539+0239. Displayed are selected regions around 
the Balmer lines, the higher Paschen series, Mg ii /1448 1 A, the Mg i b and the near-IR O i triplets. 



TABLE 1 
Results of the spectroscopic and kinematic analysis of J1539+0239. 



Vimagr 


15.72 ±0.02 


£(B-V)(mag)l' 


0.04 ± 0.03 


/i,jCOS(5(masyr"') 


-10.6 ± 1.6 


/J5(masyr"') 


-10.0 ±2.3 


/(deg) 


8.9836 


fo(deg) 


+42.9515 


/i/cos/)(masyr"') 


-14.3 ±2.1 


,Ui(masyr"') 


+2.8 ± 1.9 


TeffCK) 


7700 ± 250 


logg(cgs) 


3.00 ±0.15 


[M/H] 


-2.0a 


[a/Fe] 


+0.4 


MIMa 


0.68 ± 0.05 


rf(kpc) 


12.0 ± 2.3 


v„d(kms"l) 


-372.6 ± 5.8 


i'GRF(kms-') 


694 +300 


Vesc(kms"') 


519 







'' The visual magnitude has been derived from SDSS magnitudes 
following J ordi et al., ( .2006.) 

*" The interstellar colour excess E(B-V) has been determined by 
comparing the observed colours to synthetic ones from the model 
spectral energy distribution. 



0.4 dex with respect to iron, which is typical for the halo pop- 
ulation. We conclude that the star is a horizontal branch star 
of Population II. All results are summarised in Table [T] 

Before proceeding to the further discussion it may be in- 
structive to take a closer look on the spectrum synthesis in 
NLTE and LTE, as this has not been done for Population II 
blue horizontal branch-stars so far. A few comparisons of 
NLTE and LTE profiles are therefore shown in Fig. |3] The 
combination of a higher T^s than typically found for Pop- 
ulation II stars and the diminished line blocking because of 
the low metal content results in a hardened radiation field 
which, along with reduced thermalizing effects because of 
smaller collision rates in the low-density atmosphere, leads to 
pronounced NLTE effects on many diagnostic lines. NLTE 
strengthening is found for the Doppler core of Ha - in 
line wi th the behaviour in cool stars (e.g. iPrzvbilla & Butle j 
l2004bl) -, while the inner line wings are weakened. The 
higher Balmer and Paschen lines show much lower deviations 
from LTE. Our calculations predict the majority of the weak 



lines to be described well by the assumption of LTE. On the 
other hand, many of the stronger metal lines, as e.g. ofCi, Oi, 
Mg i/ii or Fe ii - the diagnostic lines at the spectroscopic reso- 
lution achieved within the SDSS - show pronounced NLTE 
strengthening. In order to reproduce the NLTE equivalent 
widths of these particular lines, abundance corrections need 
to be applied of about 0. 1 dex (Mg ii), 0.2 dex (Mg i, Fe ii), 
0.5 dex (Oi) and L9dex (Ci) in LTE. Note, however, that 
these line are close to saturation. LTE computations with 
increased abundances can therefore not reproduce the NLTE 
line depths at all, instead stronger line wings develop. Abun- 
dance studies based on equivalent widths may therefore be 
misleading, as such differences remain unnoticed. An inves- 
tigation at high spectral resolution would therefore be worth- 
while in order to facilitate the NLTE effects to be studied in 
detail. 

4. DISTANCE, KINEMATICS AND ERRORS 

Using the mass, effective temperature, gravity and 
extinction-corrected apparent magnitude we derive the dis- 
tance following Ramspeck et al. (2001) using the fluxes 
from the final m odel spectrum in analogy to previous work 
dTillich et alj|2009.) . The distance error is dominated by the 
gravity error. 

Applying the Galactic potential of lAllen & SantillanI 
(1991) we calculated orbits and reconstructed the path of 
the star throug h the G alactic halo with the program of 
■Odenkirchen & Broschd 0992 ). The distance of the Galactic 
center from the Sun was adopted to be S.Okpc and the Sun's 
moti on with respect to the loc al standard of rest was taken 
from lDehnen & Biniievl (119981) . As the RV is well known the 
eiTor of the space motion is dominated by that of the distance, 
ruled by the gravity error, and those of the proper motion 
components. Varying these three quantities within their re- 
spective errors we applied a Monte-Carlo procedure to derive 
the median GRF velocity vgrf at the present location and the 
velocity distribution (see Fig. |4|i and compare with the local 
escape velocity Vesc as calculated from the Galactic potential 
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Fig. 3. — Comparison of NLTE (full lines) and LTE (dotted) line profiles 
from calculations for our adopted model parameters for J1539+0239. Dis- 
played are examples of strong lines that can be used for diagnostics at the 
SDSS spectral resolution. 



of lAllen & SantilladdTool . 

A median GRF velocity of 694 km s' (with the veloc- 
ity distribution ranging from -221 to +300 km s"' around 
this value) makes J1539+0239 the fastest halo star known, 
superseding CS 221 83-0014 (at VGRF = 635+127kms-', 
iSakamoto et al.l 120031) . This value is above the local es- 
cape velocity of Vesc ~ 5 19 km s' in the potential of 
I Allen & SantillanI ( 11991b . This three component potential 
consists of a central bulge and a disk, which have a combined 
mass of Mbuige+disk ~ lO" Mq. The halo out to 100 kpc is as- 
sumed to have a total mass of Mhaio ~ 8 x 10" Mq. This total 
halo mass is insufficient to keep the star bound to the Galaxy. 

J1539-I-0239 is located in the northern Galactic hemisphere 
(/ ~ 9.0, b X! 4 2.95), in the direction clo se to the Sagittarius 
stream (Beloku rov et al.l 120061: iFellhaue r et al. 2006) and in 
particular close to the globul ar clusters NGC 590 4 and Palo- 
mar 5, with its tidal tail (Ode nkirchen et al]|2003h . However, 
it is unlikely that J1539-(-0239 is related to either of those be- 
cause of its distinct space motion and its position in the fore- 
ground. Its position and kinematics also rule out a scenario of 
a recent acceleration by an interaction with the super-massive 
black hole in the Galactic center , the p rincipal mechanism for 
generating hyper-velocity stars (lHillslll988i) . The star is cur- 
rently approaching the Galactic disk and will pass the Galac- 
tic center at a minimum distance of about 8 kpc in the future. 
See Fig. |5]for a visualization of the orbit of J1539+0239 in 
the Galactic halo (upper panel) and a for a magnification of 
the Galactic disk region of the trajectory (lower panel). 

Hence, in order to keep the trajectory of J 1539+0239 bound 
to the Galaxy, the dark matter halo mass needs to be ad- 
justed. We carried out numerical experiments increasing the 
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Fig. 4. — Galactic restframe velocity distribution for J1539+0239 derived 
from a Monte Carlo procedure with a depth of 10,000. The grey shaded 
area indicates the fraction of trajectories that would be bound for Mhaio = 
1.0 X lO'^ Mo, the most likely halo mass of lXue et alJ (20081). 



halo density by a constant factor Finally we found a bound 
trajectory for a minimum mass of MjJ™ - 1 .7 x lO'^ M©. The 
last pericenter passage occured at a distance of ~7.7kpc and 
the apocenter distance of the star's trajectory is located far out 
in the halo, at ~250 kpc in this case. If we take into account 
the full velocity distribution (see Fig. |4]i of the star we even 
can derive solutions for the extrema, which correspond to the 
absolute errors, giving Mhaio ~ l-^^fi x 10'^ Mq. 

Whether the star is bound to the Galaxy highly depends on 
the Galactic potential adopted, in partic ular on the mass o f 
the dark matter halo, as pointed out by lAbadi et alJ (l2009h . 
Our Mhaio i s similar to values found in several recent stud- 
ies. Wilkin son & EvansI (11999) used 27 satellite galaxies and 
globular clusters, by assuming that they are bound and de- 
rived a total Galactic halo mass of Mhaio ~ 1.9;^j-^ x 10'^ Mq. 
This value matches our derivation but has a larger uncertainty. 
ISakamoto et alj (l2003h used 11 satellite galaxies, 137 glob- 
ular clusters and 413 field horizontal branch stars to derive 
a total Galactic mass. The exclusion of Leo I from their 
sample would lower the total Galactic mass from Mtotai ~ 
2.5+°-^ X lO'^Mo to a value of M,o,ai ~ 1.8+^-^ x lO'^Mg, 
both again in excellent agreement with our finding. 

On the other hand, less than 4% of the trajectories resulting 
from our MC simulations would be bound for the most likely 
mass of Xue et al., Mhaio = 1 .0 x 10'^ M© (grey shaded area in 
Fig-Si- Their RV study of ~2400 blue horizontal-branch stars 
- which includes J1539+0239 but lacks any proper motion 
measurements - therefore likely underestimate the Galactic 
halo mass. 

5. SUMMARY AND CONCLUSIONS 

We reported the quantitative spectral analysis of a high- 
velocity s t ar fro m the sample of faint blue halo stars of 
IXue et alj (|2008[) . J1539+0239 was confirmed to be a Pop- 
ulation II blue horizontal branch star with a low metallicity 
of [Fe/H] - -2.0 and the characteristic enhancement of a- 
elements. Hereby, we performed a NLTE analysis of a halo 
BHB star for the first time. While the majority of the weak 
lines were confirmed to be formed close to LTE conditions, 
many of the stronger metal lines - which are important diag- 
nostics at the spectral resolution achieved within the SDSS - 
show pronounced NLTE strengthening, with the diff^erences 
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Fig. 5. — Upper panel: Trajectories for the metal-poor horizontal branch 
star J1539+0239, relatively to the Galactic disk (light blue). Applying a 
standard potential (Allen & Santillan 1991) the trajectory is unbound (black, 
t a; ±0.5 Gyr), while increasing the halo density we found a bound trajectory 
(red, -3 Gyr < ? < +2 Gyr). For reference, the Magellanic Clouds (blue dots) 
and the current position of J1539+()239 (black triangle) are marked. Lower 
panel: schematic 2D-visualization of the the central region. Details of the 
trajectory in the immediate past and around the previous pericenter passage 
of J 1539+0239 (bound solution) are shown in theX-y, Y-Z and X-Z planes. 
The positions of the Galactic center (black dot) and of the Sun are indicated. 



between the derived LTE and NLTE abundances amounting 
to 0. 1 dex to 0.5 dex typically. In addition to information on 
the chemical composition, the radial velocity, proper motion 
and spectroscopic distance were derived and a detailed kine- 
matical analysis was performed. 



Carrying out kinem atical numerical experime nts using the 
Galactic potential of lAllen & SantillanI ( 1199 Ih in order to 
obtain an orbit of J1539+0239 gravitationally bound to the 
Milky Way, we found that the mass of the dark halo has to be 
at least Mhaio ~ l-"^-!! ^ 10'^ M0(absolute uncertainties from 
extrema in MC error propagation) . This mass limi t is in good 
agreement with several previous studies (Wilkins on & Evans! 
[19991 ISakamoto et al]|2003l: lAbadi et alj2009.). Howev er, the 
significantly lower most Ukely mass value of lXue et"al] (l2008h 
is consistent with our analysis only at a level of 4%, i.e. it 
likely underestimates the Galactic dark halo mass. 

We conclude, that if the kinematics of a halo star is extraor- 
dinary enough, and the errors within the analysis are small, 
even one star alone can provide a significant lower limit to the 
dark matter halo mass, and to the total mass of the Milky Way 
(halo + bulge + disk), here Mtotai > l-8!f ] x 10'^ Mq. The de- 
termining factor is that full kinematic information is available, 
as it will become routine in the era of the Gaia space mission, 
at much higher precision. 
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